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We review the exact results for microscopic Dirac operator spectra based on either Random Matrix Theory, or, 
equivalently, chiral Lagrangians. Implications for lattice calculations are discussed. 



1. Introduction 

Exact statements can be made about the dis- 
tributions of the smallest Dirac operator eigen- 
values in gauge theories with spontaneous chiral 
symmetry breaking. This extends previous ex- 
act results for fermion zero modes in topologi- 
cally non-trivial gauge field backgrounds to an 
infinite series of non-zero modes. In this short re- 
view we focus exclusively on these exact results 
for the smallest Dirac operator eigenvalues in four 
dimensions, and discuss why they are so impor- 
tant for lattice gauge theory. 

The one single assumption that has to be made 
is that chiral symmetry is spontaneously broken. 
Actually, this statement needs a small qualifica- 
tion because the group of chiral symmetries de- 
pends on the number of fermions, and on the 
representation carried by them. We can be com- 
pletely general and consider Nf (Dirac) fermions 
in an arbitrary representation r of the gauge 
group Q with vector like couplings. As one could 
suspect from Wigner's classification of represen- 
tations, there are just three symmetry breaking 
classes to consider. It all depends on whether 
the representation is complex, real or pseudo-real. 
The case of complex representations is not only 
the physically most relevant (QCD belongs to this 
class), it is also the most simple. In the two 
other cases the initial symmetry is larger than one 
would naively have expected, namely SU(2A^/) 
rather than just SUL(iV/)xSUi?(iV/). The ex- 
pected patterns of chiral symmetry breaking are 
as follows, depending on the representation r: 
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• Complex: SV L{Nf ) xS\J R{Nf) S\J{Nf). 

• Pseudo-real: SU(2iV/) Sp(2iV/). 

• Real: SU(2iV/) -> S0{2Nf). 

By Goldstone's theorem, the cosets of the 
above symmetry breaking patterns determine the 
low-energy properties of these theories. Note that 
precise details about the gauge group Q do not en- 
ter at all. We are here getting a first bite of the 
universality that turns out to govern the lowest- 
lying Dirac operator eigenvalues. 

As explained by Verbaarschot there ex- 
ists a remarkable relation between the three sym- 
metry breaking classes above and the classifica- 
tion of chiral Random Matrix Theories ensem- 
bles. They are labelled by the so-called Dyson 
indices (3 as follows: The coset SU(7V/) of com- 
plex fermion representations corresponds to the 
chiral Unitary Ensemble chUE (/? — 2), the coset 
SU(2A^/)/Sp(2A^/) of pseudo-real representations 
corresponds to the chiral Orthogonal Ensemble 
chOE (/3 = 1), while the coset SV{2Nf)/SO{2Nf) 
of real representations corresponds to the chiral 
Symplectic Ensemble chSE of 5 = 4. This ex- 
hausts the Dyson classification 

These symmetry breaking patterns concern 
continuum fermions. Staggered fermions, away 
from the continuum, have unusual patterns of chi- 
ral symmetry breaking: all real and pseudo-real 
representations are precisely swapped. This has 
been known for a while for special cases (see, e.g. 
ref. 1^), and has very recently been shown to 

^The labelling in terms of the integers /3 is related to a 
certain number occuring in the three different Random 
Matrix Theories; there is thus no "missing" class of /3 = 3! 
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hold in general Another peculiarity of stag- 
gered fermions is their apparent insensitivity to 
gauge field topology, which we shall return to in 
detail below. 

The relationship to Random Matrix Theory 
has caused some confusion, although it can be 
stated very clearly. To get there, the theory 
first has to be formulated for fixed topological 
charge v. The order parameter for spontaneous 
chiral symmetry breaking is formally defined by 
E = lim^— >o liniy^cxD KV'V')!- To zoom in on the 
smallest Dirac operator eigenvalues, and to see 
how the chiral condensate is formed, one takes 
the chiral limit in an unorthodox way First, 
the euclidean four-volume V is kept finite, and 
it is convenient to work at fixed ultraviolet cut- 
off A. We are thus temporarily working with 
bare quantities which will have to be renormal- 
ized eventually. Both the finite volume and UV 
cut-off are ideal set-ups from the point of view of 
lattice gauge theory. 

Next, the quark masses are chosen so that 
the pseudo-Goldstone bosons ("pions") are too 
light to fit inside the box, i.e. TOtt ^ ^/L^ 
where L is the linear extent of the 4-volume V . 
Now V is sent to infinity while the combina- 
tion rriiTiV {rrii being quark masses) is kept fi- 
nite - it is always possible to satisfy these two 
constraints by taking mi small enough. Because 
L is taken much larger than \/Kqcd, the eu- 
clidean partition function is dominated by the 
pions. Higher- mass states do contribute to the 
partition function, but as the volume is sent to in- 
finity their contribution will eventually, for large 
enough volume, be exponentially suppressed, of 
order exp[— ML], with ML ^ 1. In contrast, 
the pions do not yield exponentially suppressed 
contributions: their masses can be tuned to zero, 
and we precisely require ttIt^L <^ 1 throughout. 
It is in this sense that the results which will be 
reviewed below are exact: given any required de- 
gree of accuracy, this accuracy can be achieved by 
simply tuning V and m^. The exact results hold 
in the limit. 

Leutwyler and Smilga proposed to combine the 
above ingredients oi V ^ 1/ and fixed topol- 
ogy Because the euclidean partition function 
is dominated by the pions, it follows from the 



coset of chiral symmetry breaking. It is a non- 
renormalizable chiral Lagrangian with, in princi- 
ple, an infinite number of terms. However, be- 
cause of the peculiar finite-volume regime cho- 
sen, only the zero-momentum modes need to be 
considered (again an approximation that can be 
made as accurate as we wish by tuning rrii and 
V) 1^. Let us focus mainly on the symmetry 
breaking class relevant for QCD. Then the effec- 
tive partition function is dominated by one single 
term in the effective Lagrangian: 



dU exp 



SU{Nf) 



VEReTiie'^ /^fMU^)] (1) 



One crucial point here is the dependence on the 
vacuum angle 9. Because of the anomaly, the par- 
tition function does not depend on the quark mass 
matrix M and 9 separately. A chiral rotation (in- 
dependent phase rotations of the left and right 
handed fields) gives a linear shift in 9, and since 
the partition function is left unchanged by such 
a change of integration variables, it can only de- 
pend on the invariant combination exp[i9/Nf]M.. 
We normally do not consider QCD with a non- 
zero vacuum angle, and here it is also just intro- 
duced at a preliminary step, as a source of topo- 
logical charge i^, and as a means for computing 
the partition function in sectors of fixed v jij : 

Z^= [ dU (detC/)''exp[yEReTr(Alt/^)] (2) 

Ju{Nf) 

The projection onto fixed topological charge has 
been used as the additional U(l)-factor that ex- 
tends the integration over the zero-momentum 
modes to U(iV/). Note that Z^, does not depend 
on mi,V and E independently, but only on the 
combination 



(3) 



We will see that this finite-volume scaling is re- 
flected in the Dirac operator spectrum. 

There are two distinct advantages for consider- 
ing the theory at fixed topological charge v. The 
first is that the zero-dimensional group integral 
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can be done analytically for all Nf and i' 



Nf 



where /„(x) is a modified Bessel function, and the 
matrix in the numerator is of size Nf x Nf. The 
second is the big surprise: at fixed z/ there is a 
connection to Random Matrix Theory. 

2. The relation to Random Matrix Theory 

As observed by Shuryak and Verbaarschot 0, 
the group integral can be rewritten as a Ran- 
dom Matrix Theory partition function which has 
an uncanny resemblance to the original QCD 
path integral. Consider 



Nf 

Z^ = y^W^J]^ det(i7\f+m/)exp 



where 



M = 





w 



-^Try(M^) 



(5) 







(6) 



Here the integral is over complex matrices W of 
rectangular size N x {N + v). The Wi's are di- 
mensionless numbers, and the potential V{M'^) is 
unspecificied at this point. The matrix M anti- 
commutes with diagjljv, — ^n+iAi ^'^'^ ^ ^ con- 
sequence the eigenvalues of M occur in pairs ±A 
whenever A ^ 0. Because of the rectangular na- 
ture of W , the matrix M also has precisely v zero 
modes. The intuitive idea is the analogy with the 
determinant of the Dirac operator for complex 
representation fermions: the matrix M has v zero 
modes, it is chiral, and in (^) one integrates over 
complex matrices. 

But one is not simply trying to replace the 
path integral over gauge potentials (x) by zero- 
dimensional matrices. Instead, the precise rela- 
tionship is as follows. Take a "microscopic limit" 
of eq. (||) in which jli = mi{2N)p{0) is kept fixed 
as iV ^ oo. In that limit, and up to an irrelevant 
(/ii-independent) normalization. 



(7) 
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This was first demonstrated in ref. for the 
case of a Gaussian potential in the (3 = 2 chiral 
ensemble. It follows from a series of universality 
theorems that the identity holds for any generic 
choice of V{M) ||^, and that a fine tuning is re- 
quired to reach other universality classes that are 
of no obvious relevance here. Basically, the do- 
main of universality is determined by the condi- 
tion that p(0) ^ 0, as could have been expected 
from the manner in which the microscopic limit 
is taken. Identities similar to exist for the two 
other classes of chiral symmetry breaking ||l^ . As 
outlined above, the chiral Random Matrix Theory 
ensembles are different, but there are analogous 
universality proofs for those cases jll] . 

Having two partition functions coincide as in 
eq. (|^) may not seem to give much informa- 
tion. But one should rather view Z^[{fj,i}] as 
a generating function for the chiral condensate. 
Then the statement (|7|) is highly non-trivial since 
it means that in the infinite-volume limit where 
fj,i = rriiTiV is kept fixed we can just as well 
compute the chiral condensate from Random Ma- 
trix Theory. This observation is also the basis 
for finally understanding why it is even possible 
to compute the microscopic Dirac operator spec- 
trum from Random Matrix Theory. More about 
this below. 

3. The microscopic Dirac operator spec- 
trum 

Although the two partition functions (||) and 
(|^) coincide in the microscopic limit, it is far from 
obvious that the eigenvalues A of the matrix M 
should be related to the eigenvalues of the Dirac 
operator A. The eigenvalue spectrum of the Ran- 
dom Matrix Theory is very different from that of 
the Dirac operator: for a Gaussian potential the 
eigenvalue density of (^) is the famous Wigner 
semi-circle, while the spectrum of the Dirac op- 
erator is expected to approach p{X) ^ \^ near 
the UV cut-off. One may search for universal- 
ity that makes such differences irrelevant. Here 
it means going to a scale near the origin where 
both spectra, on that scale, are trivially identical 
(namely constant, near A ~ 0). This is possible, 
because by the Banks-Casher relation E = 7rp(0) 
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and the assumption of spontaneous chiral symme- 
try breaking we know that p(0) ^ 0. Similarly, 
we must require of the potential in that it 
leads to p(0) ^ 0. Then one can define micro- 
scopically rescaled variables, which for the Dirac 
operator eigenvalues are C = AI^S (cf. eq. (||)), 
and for Random Matrix Theory C = A27V7rp(0). 
A microscopic spectral density (and similarly for 
higher spectral correlators) of the Dirac operator 
is analogously Ps(C) = p(C/(^S))/^i which has 
a finite well-defined limit as — > cxo 0| . The mi- 
croscopic density in the Random Matrix Theory 
context is defined analogously, and it can be com- 
puted by various technqiues. The first analytical 
expression, for TV/ massless flavors in the (3 = 2 
universality class, was obtained in ref. p^ ]: 

C 2 

Ps{0 = -^[Jn f+iy{C) ~JNf+u-l{C)JNf+u+l{0] (8) 

which, as a first check, was found to reproduce 
the Leytwyler-Smilga spectral sum rule for e.g. 
(Snl/Cn) [ P^IH]- Corresponding expressions ex- 
ist for P = 1,4 |ll|]. 

Quark masses can be included without diffi- 
culty by performing a "double-microscopic limit" 
in which both rrn and eigenvalues A are blown up 
as in (^). The precise analytical form of the spec- 
tral correlators can then be worked out explicitly 

There is a very compact formulation of all these 
results. A central object is the so-called kernel 
K{(a,Cb) from which all spectral correlation func- 
tions can be computed: 



PsiCl 



detX(Ca,Cfc) 

ab 



(9) 



A "Master Formula" for this kernel is |]T2| 

Nf 



z^''^^'\{^,},^C,^C)/zi''^\{^,}) . (10) 

This formula encompasses all cases for = 2, any 
number Nj of fermions with masses /i/, and any 
v. Moreover, although it is derived within Ran- 
dom Matrix Theory, we have used the identity (0) 
to express the r.h.s. entirely in terms of field the- 
ory partition functions. There are closely related 
formulas for the /3 = 4 case p2| . 



One can also compute an infinite series of in- 
dividual eigenvalue distributions, beginning with 
the smallest. They too can be expressed in terms 
of the effective partition functions [Q . The ana- 
lytical formula is known in all generality: it gives 
the /c'th eigenvalue distribution for all three uni- 
versality classes, for any number of fermions TV/, 
and in an arbitrary sector of topological charge 
To give a simple example, the quenched dis- 
tribution of the first eigenvalue in a sector of ar- 
bitrary u reads (i, j = 1, . . . , i/) ll^ 



^min(C) = ^e-«'/4det[/2+.-,(C)] 



(11) 



for the (3 = 2 universality class. 

One can perform the sum over topological 
charge explicitly for all spectral correlation func- 
tions, including the microscopic spectral density 
itself Q. The result can for Nf > he writ- 
ten in terms of the effective field theory partition 
functions at vacuum angle = n when k is odd: 

ft(Ci,..^.a) = (-i)'*'''"n(c|-c?)' X 

n(Kin«?.M;)) "'"':j!;;':!:«'"- (i2) 



4. Derivation from chiral Lagrangians 

All microscopic spectral correlation functions, 
and all probability distributions of individual 
eigenvalues, can thus be expressed in terms of 
effective field theory partition functions. This 
puts universality of these quantities on a very 
simple footing. It also gives a strong hint that 
in fact Random Matrix Theory is not needed 
at all. These results can indeed be derived en- 
tirely within the chiral Lagrangian framework. It 
hinges on the crucial identity (^, and an exten- 
sion of it which we will discuss next. 

One way to compute the spectral density of the 
Random Matrix Theory (^) is through the resol- 
vent 



Ri,{m.o) = ( Tr 



1 



M - 



m, 



(13) 



*The only exception is tlie /3 = 1 class where, for technical 
reasons, f must be odd in the general formula. 
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where rriv is an external parameter. In the mi- 
croscopic hmit this is, in view of eq. (Q), nothing 
but the partially quenched chiral condensate. To 
be precise, consider the supersymmetric method 
for computing this quantity from the chiral La- 
grangian |2^: One adds a quark of mass my + j, 
and a bosonic quark of mass mv Then the par- 
tially quenched chiral condensate is dj lnZ^|^-^Q. 
In a similar way, one can modify the Random 
Matrix Theory and multiply the integrand by 
det[A/ -I- TO„ -I- j']/det[M -I- m„]. Taking the mi- 
croscopic limit, this leads to a supersymmetric 
generalization of the identity (0), where the chi- 
ral Lagrangian is the supersymmetric one asso- 
ciated with partial quenching. Naively this La- 
grangian would be based on the Goldstone su- 
permanifold U{Nf + 1|1), but in fact the proper 
choice of integration domain is a more subtle 
issue, as discussed at length in ||2^,|2^. This 
way of proceeding is only in order to establish 
why the Random Matrix Theory approach really 
does give exact answers for the microscopic Dirac 
operator spectrum. If one is not interested in 
this, one can proceed completely within the chi- 
ral Lagrangian framework, and simply compute 
the partially quenched chiral condensate directly. 
As a simple (3 = 2 example, consider the fully 
quenched case where, with fj,^ = fj,, one finds [pT| : 

M^=M(/.(A*)i^.(M)+/.+i(M)^.-i(M))+-(14) 

The microscopic spectral density of the Dirac op- 
erator is given by 

Ps(C)-^DiscI],(/i)|^^,;^ (15) 

which precisely agrees with 

Similar expressions exist for higher spectral k- 
point functions. See ref. for the (3 = 2 case. 
There is again exact agreement with the Random 
Matrix Theory results. 

5. Replicas 

It is interesting that one can get quite far by 
an alternative formulation of partial quenching 

®The microscopic spectral density by definition does not 
include the zero modes. 



based on the replica method. The fermion deter- 
minant is then removed by analytic continuation 
in the number of flavors, so that eventually one 
can take the limit 7V„ 0. Here iV„ denotes a 
number of additional unphysical flavors added to 
the theory with Nf physical quarks. 

The first observation is that now the identity 
as it stands is all that is needed in order to 
show equivalence to the Random Matrix Theory 
approach. The integration in eq. (Q) is over 
U(iV/ -I- iV„), and the limit ^ is taken 
after differentiating w.r.t. fiy. Small-mass and 
large-mass expansions can be performed in this 
way |2^,^. As an example, for the quenched 
chiral condensate the replica method yields pre- 
cisely small-mass and large-mass series expan- 
sions of the expression whose inversion in 
turns yields the same microscopic spectral den- 
sity as was originally derived from Random Ma- 
trix Theory. It is a highly non-trivial check on 
the consistency of the whole formalism that these 
results all coincide. 

6. Lattice results 

This review has not proceeded in historical or- 
der. The derivations of the microscopic Dirac op- 
erator spectrum based directly on the partially 
quenched chiral Lagrangian came after the first 
lattice study already had demonstrated the 
viability of the Random Matrix Theory approach. 
Since then there have been numerous lattice com- 
parisons with the exact analytical expressions, 
for different gauge groups, different representa- 
tion quarks in both quenched and full theories 
|27| , |2^ ,p[| . It is impossible to do justice to all that 
work here, but a few points can be highlighted. 

A quite systematic study of the smallest eigen- 
value distributions in all three universality classes 
and for zero and non-zero topological charge 
showed perfect agreement with the analytical pre- 
dictions ||2^. The distribututions of a whole se- 
quence of smallest eigenvalues have been com- 
pared with the analytical formulas in ref. 
One sees clearly how the microscopic spectral 
density is built up of the individual distributions. 

Also higher-point spectral correlators have 
been checked again with spectacularly 
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Figure 1. The first Monte Carlo measurement of 
the microscopic spectral density of the Dirac op- 
erator for gauge group SU(2) and quenched 
staggered fermions. The dashed curve is the an- 
alytical prediction. 



good agreement with the analytical predictions. 
Extensions to the regime of chiral perturbation 
theory have also been made [32[ . 

Lattice calculations with staggered fermions 
compare only with the analytical expressions for 
the sector with i' — Q. This makes sense, since 
it has been known for long ||3^] that staggered 
fermions do not have exact zero modes at usual 
gauge couplings. Even if one selects distinct topo- 
logically non-trivial gauge field sectors, the mi- 
croscopic Dirac spectrum seems completely unaf- 
fected, and always agrees with the v = predic- 
tions |34[ . There is a much more direct way of 
understanding this. The symmetry breaking pat- 
terns of staggered fermions away from the con- 
tinuum always have an additional U(l) factor^ 
which for the zero momentum modes is com- 
pletely equivalent to a projection on the 1^ = 
sector js^,^. So we should only compare with 
what can be viewed as v — predictions of the 
theory without the extra U(l) factor. If one sees 
sensitivity to topology [p6| , then continuum flavor 



^They are V(N)xV{N) - 
Sp(27V) [real], and U{2Af) 



symmetries are beginning to be recovered. 

Instead of looking directly at the smallest Dirac 
operator eigenvalues, one can also consider de- 
rived quantities such as the (mass-dependent) chi- 
ral condensate p7|-pl|| 



= 2/i 



dc 



or higher chiral susceptibilities, such as 



dC 



.2\2 



2v 
~7i 



(16) 



(17) 



U(7V) [complex], U(2Af) 
• SO(2Af) [pseudo-real] 



/O '(C^+M^J^ 

In the quenched and partially quenched cases 
there is a one-to-one correspondence between 
S(/x)i, and the microscopic spectral density 
p^s'^{Q. But: ^u{p) has a short-distance singu- 
larity that is not taken into account in the above 
description. At fixed UV cut-ofF A this is not a 
problem when one takes the infinite- volume limit, 
as one should in order to compare with ([l^). 
Corrections are of the form AmA^ + _Bm'^ In A, 
where A and B are constants. This looks scary, 
but we are taking the limit in which /x = mTiV 
is kept fixed as F — > oo. So these UV cor- 
rections actually vanish; they are suppressed by 
inverse factors of the volume, and really read 
{A/^){^l/V)A^ + {B/T,^){fxyV^)lnA. If the vol- 
ume V is not large enough these terms can be 
annoying, and then one should subtract them. 

There are cases where very, very few Dirac op- 
erator eigenvalues play any role in building up 
the chiral condensate ([ij) in the interesting mass 
range. As an example, let us compare the exact 
result dl^ ) with what one gets by keeping only 
one single Dirac operator eigenvalue in the inte- 
gral of eq. ( p^ for v — 0. We then simply replace 
p^s^\0 ill the integrand of ( p^ by the probabil- 
ity distribution of just one single Dirac operator 
eigenvalue, the smallest, P^[^iC) = (C/2)e~''^/2. 
This may seem like an absurd truncation, and 
of course an unnecessary one, since we know the 
full analytical form of p^s'\(). But the compar- 
ison between this approximation, and the exact 
answer is shown in figure 2. 

The two curves are basically indistinguishable 
(up to fi ^ 0.1 the difference is less than 5%, de- 
creasing fast as /i is lowered). If one sees good 
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Figure 2. The exact quenched chiral condensate 
(here m stands for /i) for 1^ = 0, compared with 
keeping only one single eigenvalue in eq. (16). 



scaling of the i/ — quenched chiral condensate 
in such a mass range, then one is probing, and 
having good statistics of, only one single Dirac 
operator eigenvalue. Of course, we could just as 
well have included all the other eigenvalues, and 
then the relation would be exact. For extracting 
S in this way it seems more logical to measure the 
distributions of the Dirac operator eigenvalues di- 
rectly. When ly ^ the eigenvalues are shifted 
further away from the origin, the integral (^) is 
not as infrared sensitive, and the brutal approxi- 
mation of keeping only the first eigenvalue is not 
nearly as good. The error is easily assessed if we 
rewrite, for v ^ p8|. 



5]^(m) 




(18) 



which shows that to leading order one is simply 
trying to measure the first Leutwyler-Smilga sum 
rule 0, (Enl/C«) = l/(4t^)- It seems advanta- 
geous to instead measure the first few eigenvalue 
distributions, or, for lower statistics data, at least 
individual averages. That way one avoids con- 
taminating the result with the larger eigenvalues 
which at the limited finite volume will have distri- 
butions that are incompatible with the analytical 
formulae. 



7. Beyond Random Matrix Theory 

We have emphasized that the results reviewed 
here are exact in the sense that they can be made 
as accurate as we wish by tuning the quark masses 
rrii and the volume V. It immediately suggests 
that it is possible to compute corrections to this 
scaling, which would improve convergence and 
bring in more terms from the chiral Lagrangian. 
This is indeed the case. A computational frame- 
work was laid out by Gasser and Leutwyler, who 
called it the e-expansion 

The idea is simple. Instead of keeping only the 
zero-momentum mode in the chiral Lagrangian, 
one includes the modes of non-zero momentum 
in a perturbative manner. Thus, one starts with 
the full chiral Lagrangian and then separates out 
the zero-momentum modes of U{x): 



Uix) 



(19) 



where u G SU(A^/) is a space-time independent 
collective field, and ^(x) contains only modes of 
non-zero momentum. To leading order, where one 
simply ignores S,(x), this gives the effective par- 
tition function (|l|) after identifying U = v?. To 
next order the kinetic energy term contributes, 
and one resorts to the usual loop expansion in the 
fluctuation field ^(x). The modes of zero momen- 
tum are still treated exactly. So this expansion 
is a combination of the exact, non-perturbative, 
leading-order result for y oo, and a pertur- 
bative expansion where each term should be sup- 
pressed by powers of 1/L = 1/]/^/^. In 4 dimen- 
sions the expansion is even better behaved, as the 
first correction goes as 1/i^, rather than just 1/L 
1^,^. Leading terms are what we alternatively 
could obtain from Random Matrix Theory. The 
corrections take us beyond. 

The one-loop correction to the chiral conden- 
sate computed in this way can be re-absorbed into 
a volume-dependent Ee//(V^) fl, 



= 1 



N]-l 



1 PiiU/L) 



N 



f 



p2 12 



(20) 



where /3i depends on the geometry ||4^,|45[|. This 

^The quantity e counts orders in 1/L, wliere L is the linear 
extent of the volume. There is no relation to the expansion 
around critical dimensions. 
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carries directly over to the Dirac operator eigen- 
values, on account of the relation ([l5|). From the 
small 1/L^-correction one can thus, with sufh- 
cient statistics, also measure the pion decay con- 
stant F from the eigenvalue distributions. 

The quenched e-expansion meets the usual dif- 
ficulties of the quenched theory. We outline 
the problem here using the replica version of 
quenched chiral perturbation theory 1 43 1 because 
it is simpler to explain. Then ( |l9| ) still defines the 
split between modes of zero and non-zero momen- 
tum, but now U{x) GU(iV„). Let S(a;) = Trf (x) 
denote the flavor singlet fluctuations. The lowest- 
order chiral Lagrangian in a sector of fixed topol- 
ogy reads [p5| 



Tr 



(21) 



and the propagator of S(a;) is thus modified due 
to the mQ-term (the precise form is given in 
p3|]), and one again takes the limit A^^ ^ 
after having differentiated with respect to the 
sources. A new scale is introduced by the mo- 
term, and in order to proceed perturbatively one 
must assume that rn^/ {^ttF'^) is small. At one- 
loop level the quenched analog of the effective 
volume-dependent Se//(y) is [|5| 



E 3F2 



(22) 



In the full theory logarithms are typically multi- 
plied by , and hence in this regime suppressed 
by at least Here it survives at fixed mo, 

indicating difficulties with the whole expansion. 
It is not a quenched chiral logarithm, but, in this 
regime, a quenched finite-volume logarithm that 
signals the limitations of the expansion. 

The e-expansion has also been worked out for 
correlation functions of the full theory 
and very recently these results have been gener- 
alized to the quenched case [Q. By measuring 
correlation functions in this regime one can ex- 
tract the low-energy constants of QCD. 



8. Conclusions 

The microscopic tail of the Dirac operator spec- 
trum is computable by a precise comparison with 
the associated effective field theory. The result 
is a series of exact expressions for observables in 
finite-volume gauge theories. This is a highly un- 
usual situation for such strongly coupled theories, 
and one can make good use of these results. In 
particular: 

• Reliability of fermion algorithms can be 
checked by comparing with the exact expressions. 

• Gauge field topology gives distinct predic- 
tions, and is probed at the "quantum level" . 

• Physical observables derivable from the mi- 
croscopic spectral correlation functions are known 
in exact analytical forms. Simple models, 
ansatze, or unknown extrapolations for the small- 
mass finite-volume behavior can be replaced by 
the correct analytical expressions. 

• By going to unphysical regimes we can ex- 
tract physical observables. This includes unphys- 
ical volumes, and, if need be, fixed gauge field 
topology. The most extreme example was given: 
From the distribution of just one single Dirac 
eigenvalue we learn i) whether spontaneous chiral 
symmetry occurs in the infinite- volume theory, ii) 
which symmetry breaking class it belongs to, iii) 
the value of the infinite- volume chiral condensate 
E, and finally, by looking even closer, iv) the pion 
decay constant F. If one is willing to push it, in 
principle the infinite series of parameters of the 
chiral Lagrangian can be probed by the Dirac op- 
erator eigenvalues. 
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